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Abstract 

The temperature anisotropies and polarization of the cosmic microwave background radiation 
(CMB) not only serve as indispensable cosmological probes, but also provide a unique channel 
to detect relic gravitational waves (RGW) at very long wavelengths. Analytical studies of the 
anisotropies and polarization improve our understanding of various cosmic processes and help 
to separate the contribution of RGW from that of density perturbations. 

We present a detailed analytical calculation of CMB temperature anisotropies and po- 
larization Pk generated by scalar metric perturbations in synchronous gauge, parallel to our 
previous work with RGW as a generating source. This is realized primarily by an analytic 
time-integration of Boltzmann's equation, yielding the closed forms of a^. and /J^. Approxi- 
mations, such as the tight-coupling approximation for photons a prior to the recombination 
and the long wavelength limit for scalar perturbations are used. The residual gauge modes in 
scalar perturbations are analyzed and a proper joining condition of scalar perturbations at the 
radiation-matter equality is chosen, ensuring the continuity of energy perturbation. 

The resulting analytic expressions of the multipole moments of polarization af , and of tem- 
perature anisotropies af are explicit functions of the scalar perturbations, recombination time, 
recombination width, photon free streaming damping factor, baryon fraction, initial amplitude, 
primordial scalar spectral index, and the running index. These results show that a longer re- 
combination width yields higher amplitudes of polarization on large scales and more damping 
on small scales, and that a late recombination time shifts the peaks of Cf^' to larger angular 
scales. 

*caiz at email.arizona.edu 
tyzh at ustc.edu.cn 
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Calculations show that af is generated in the presence of the quadrupole a2 of temperature 
anisotropies via scattering, both having similar structures and being smaller than the total af, 
which consists of the contributions from the monopole, dipole, quadrupole, and Sachs-Wolfe 
terms as well. The origin of the two bumps in Cf'^ on large angular scales is found to be due 
to the time derivative of the monopole of temperature anisotropies. Furthermore, af together 
with af demonstrates explicitly that the peaks of Cf^ and Cf^ alternate in I— space. These 
results substantially extend earlier analytic work. 

The analytic spectra C^^' agree with the numerical ones and with those observed by WMAP 
on large scales {I < 500), but deviate considerably from the numerical results on smaller 
scales, showing the limitations of our approximate analytic calculations. Several possible 
improvements are pointed out for further studies. 

PACS number: 98.70. Vc, 98.80.-k, 98.80.Jk, 

Key words: cosmic microwave background radiation, scalar perturbations, polarizations 
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1. Introduction 



By confronting predictions of theoretical cosmological models with the data on the CMB by 
the observations, such as BOOMERANG [T], MAXIMA [2], DASI [3], WMAP [3 E El [3 [E E], 
Archeops [lOj, CBI [TT|, QUaD [12J, BICEP [13] etc, several important cosmological parameters 
of the standard Big Bang model have been directly measured or constrained. These studies have 
been instrumental for rapid progresses toward understanding of the evolution of the Universe, and 
for the advent of an epoch of precise cosmology. 

On the side of theory, these achievements have been possible through detailed computations of 
the spectra of the CMB temperature anisotropies and polarizations. Codes of numerical compu- 
tation, such as CMBFAST [14] and CAMB [15j, give the the spectra Cf^' of CMB temperature 
anisotropies and polarizations. The prominent structure of C{^^' involves various cosmological 
parameters, as it depends upon several major physical processes during the cosmic expansion, such 
as the inflation, radiation-matter equality, recombination, and the reionization as well. Analytical 
studies are still indispensable for understanding how various underlying physical effects give rise to 
the observed behavior and for theoretical interpretations of the observational data. In particular, the 
analytical spectra are helpful in revealing their explicit dependence on the cosmological parameters 
and possible degeneracies between them. So it would be desired to have the analytical C^^^' for a 
better understanding of physics of CMB. 

From computational point of view, the CMB temperature anisotropies and polarizations are 
determined by the Boltzamann's equation of the photon gas in the expanding Universe. Although a 
number of ingredients will influence this equation, two key elements are responsible for the overall 
features of Cf^', i.e., Thompson scattering during the recombination process around a redshift 
z ~ 1000 and the metric perturbations hij of Robertson-Walker spacetime entering the equation as 
the Sachs- Wolfe term [16] . Generally, there are two types of metric perturbations as the source: the 
scalar (density) perturbations [171 HE [IS] and the tensorial perturbations, i.e., RGW [23 EH 122] ■ 
Both types can be generated during early stages of the universe, such as the inflationary expansion. 
Among them, the contribution from scalar perturbations is believed to be dominant over that from 
RGW [23l[24ll25ll26l[2Zll28ll29ll30], characterized by a tensor/scalar ratio r [SB [32]. For a power 
law spectrum of the primordial fluctuations, WMAP5 data alone puts an upper limit on the ratio 
r < 0.55 (95% CL) [7j, while WMAP7 gives r < 0.49 (95% CL) for ACDM+Tensors+Running [8]. 
The recent data of LIGO S5 [33] with cross-correlation of HI and LI gives a constraint r < 0.55 for 
the flat primordial tensorial perturbations with a negligible running index [34]- For the case of RGW 
as the source, Refs. [351 [35] derived the analytical spectrum Cf^ of temperature anisotropies, and 
Ref. [37] gave all four analytical spectra Cf^, C^^^, Cf^, as well as Cf'^. Ref. [3B] incorporated 
the reionization process into calculation and obtained the reionized analytical spectra C^^' . 



3 



Ref.[35] presented a fully covariant and gauge-invariant formulation of Boltzmann's equations. 
The analytic calculation of the scalar induced Cf^ was made in Newtonian (longitudinal) gauge 
in Refs |4U1 |4l]. Ref.[42j gives an analysis of Cf^ in synchronous gauge, but the treatment 
of temperature anisotropies itself was not enough to separate contributions of monopole, dipole, 
quadrupole, and Sachs-Wolfe terms. It did not address the CMB polarizations either. Ref.[43] 
gave a unifying framework in synchronous gauge to discuss the scalar induced spectra Cj^ , Cf^, 
Cf^ and as well as the RGW induced spectra C^^' . Motivated by possible extractions of signals 
of RGW using anti-correlation of Cf^^ , attempts were made to estimate qualitatively the possible 
forms of multipoles of of the temperature anisotropies and af of the polarization at / ~ 50 [43] . 
However, the analysis was still preliminary by lacking of an explicit formula of af , since the time- 
integrations of the Boltzmann's equation as a key procedure was not carried out. Viewing these, 
in this paper, we shall perform a detailed, analytic calculation of af and af induced by the scalar 
perturbations in synchronous gauge, and present the analytical spectra Cf-^' , which will be at the 
same level of accuracy as the analytical Cf^' by RGW [361 [371 [38]. Aside several new insights 
into the physics of CMB, in particular, our resulting cross-correlation spectrum Cf^ has already 
demonstrated some inaccuracy in the preliminary analysis of Ref.[43j. Therefore, these two sets of 
analytic spectra together are more reliable in analyzing and disentangling the RGW contributions 
from observational data [44l l45l l46]. 

The synchronous gauge has been often used, in which the decomposition of generic metric 
perturbations hij into the scalar, vector, and tensorial types is straightforward. For the scalar 
metric perturbations, this gauge is also more efficient in dealing with the adiabatic and isocurvature 
initial conditions, adequate for numerical computations [14l|T5]. In comparison with the conformal 
Newtonian gauge [40], there are residual gauge freedoms in the synchronous gauge in the solution 
of scalar metric perturbations. By the restricted coordinate transformations [471 [481 [49] , general 
solutions become rather involved for modes of arbitrary wavelengths. To implement analytical 
calculations, we work in the long wavelength approximation. Besides, a joining condition of the 
perturbation modes at the equality of radiation-matter will be chosen to ensure the continuity of 
the energy perturbations, not of the pressure. 

In solving the Boltzmann's equation. One has to carry out the time-integrations for af and 
af from the RD epoch up to the present. The visibility function for the recombination process 
will appear in the integrations, and can be approximately fitted by the Gaussian type of functions 
[361 [371 [381 [50] , in order to obtain the analytical expressions of af and af. 

The organization of this paper is as follows. In Section 2, we introduce the convention of the 
decomposition of the scalar metric perturbations hij into two independent modes /i(r) and r](r) in 
the flat Robertson-Walker metric. In Section 3, the Boltzmann's equation of the CMB radiation 
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field in the Basko-Polnarev's framework is formally solved in terms of two time-integrations for 
the temperature anisotropies akir) and polarizations /3fc(r), respectively. The integrands consist 
of some combinations of the metric perturbations, the monopole ao, and the dipole ai of the 
temperature anisotropies as well. The fitting formula for the visibility function involved in the 
integrand is introduced. In Section 4, in the tight-coupling approximation, both ao and ai are 
solved in terms of the metric perturbations. In Section 5, the time-integrations are carried out, 
yielding the analytical expressions of af (r) and af (r), respectively. In Section 6, we will remove 
the gauge modes from the scalar metric perturbations for the RD and the MD eras, make a joining 
connection of the perturbations at the radiation-matter equality, and choose the proper initial 
conditions for the perturbations. In Section 7, we present the final analytical spectra Cf^^, Cf^^, 
and Cf^^, and compare them with the numerical and the observed results. Several interesting 
properties of CMB anisotropies and polarization are revealed by the analytic spectra. Section 8 
summarizes the main results and discusses possible future improvements. The Appendix provides 
the formulae that relate the multipole moments of and af to and (3^, respectively. The unit 
with h = c = = I will be used. 

2. Scalar Metric Perturbations in Synchronous Gauge 

For a spatially flat {k = 0) Robertson-Walker (RW) space-time, the metric is 

ds^ = a\T)[-dT^ + {6,j + hij)dx'dx^]. (1) 

where a(r) is the scale factor as a function of the comoving time r. The normalization of the scale 
factor is taken such that a(ro) = 2/Hq at the present time tq, where Hq is the Hubble constant. 
In our calculation, the RD and MD stages of the cosmic expansion are involved, for which the scale 
factor can be taken as the following form 
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air) = ^iT + T2Y, T>T2, (3) 







respectively, where T2 = l/{2y/l + z^q) — 0.0085 is the radiation-matter equality for Zeq — 3400 
|4], and tq = 1 — T2. In this convention, the recombination time ~ 0.0216, corresponding to 
a redshift z ~ 1100. As will be seen in Section 3, the precise value of actually depends upon 
the baryon fraction Qh- To keep our analytical calculations simple, we do not include the current 
accelerating stage, which will bring some minor modifications to the CMB spectra. The metric 
perturbations in the synchronous gauge hij in Eq.([T]) can be generally decomposed as 

h^J = lh6,, + hi + hj-j + hi . (4) 



Here /i^ is the transverse {didjhf- = 0), vector mode, and is usually neglected as it decays with 
the cosmic expansion, hjj is the transverse {dihjj = 0), tensorial mode, i.e., RGW. Its analytic 
solution and the analytical spectra Cf^' induced by hjj have been studied before SZl SHI E] . 
We consider the remaining part, which is the scalar metric perturbations, 

h^j = ^h6ij + hi , (5) 
where h = h\\s the trace part, and h}]- is the traceless and longitudinal part, satisfying 

eijkdj dihl = 0. 
It can be expressed in terms of a scalar function, 

4 = (d^d, - U,vA V. (6) 
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Thus the density perturbations hij are described by two scalar functions, and can be written as a 
Fourier integration [431 122] 



(7) 
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where hkiT) and hkiT) are the two scalar functions introduced, and 
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Pij = 4,%, Pij = (4,% - gf^ij) 

are the two corresponding polarization tensors for the density perturbations. If we write 

1 2 

hk{r) = hkij), hkij) = Qrikij), (9) 
then /ifc(r) and %(r) are the scalar modes used in Ref.[52j. If we write 

hk{r) = -Hi fc(r) + SHkir), hk{r) = -SHj^ir), (10) 

then Hk and Hi ^ are identified as the scalar modes employed in Refs. [451 112], where small letters 
hk and hi k were used. The sets {h,r]) in Eq.(l9]) is related to the set {H,Hi) in Eq.lITOl) as the 
following 

H = -2r], Hi = -{h + 6T]), (11) 

where the sub-index k has been omitted in the following when no confusion arises. 

An important property of density perturbations is that, a k mode of hij in Eq.(|5]) is rotationally 
symmetric about the k axis. Let the polar axis z be along k. The k mode of the trace part ^hSij is 
isotropic in space, and the longitudinal part h^ has only the zz component. This is also reflected 
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by its polarizations Pij and Pij given in Eq.(8) that only depend on the vector k, independent of 
any vector perpendicular to k. So the k mode of hij is independent of the azimuthal angle (p. 
In contrast, a k mode (plane wave) of GW hjj is transverse, with two components h^x = ~hyy 
and hxy = hy^. Therefore, hjj is not rotationally symmetric about the k axis, and does depend 
on the azimuthal angle 0. Due to this difference, as we shall see later, the 0-independent density 
perturbation does not induce the magnetic type of polarization of CMB, whereas the ^-dependent 
GW does. 

In order to calculate the evolution of CMB anisotropies and polarization, one needs the dynamic 
evolution of scalar perturbations /i(r) and r](r) that enters the Boltzmann's equation of photons. 
However, in the synchronous gauge, the solution of /i(r) and r](r) contain the residual gauge modes 
for both the RD and MD stages, which have to be dealt with later (in Section 6). 

3. Boltzmann's equation in RW spacetime 

The temperature field of CMB is not exactly isotropic, instead it has anisotropies, which are 
related to the metric perturbations hij via the Sachs-Wolf term. Moreover, the quadrupole com- 
ponent of the temperature anisotropies will be further induced by the linear polarizations via the 
Thomson scattering during the recombination. So the radiation field is described by the following 
column vector [431 [231 |24l [531 IM] 




^=7^-, I-Q , (12) 



where / is the intensity of radiation, and Q and U together describe the linear polarizations. The 
column n can be split into two parts 



n 




n^^^u + n^^'^ with u = I 1 , (13) 



where n^^^ is the homogeneous, isotropic and unpolarized Planck spectrum in the expanding universe 
with frequency rescaled by the scale factor u = ua^r), and n'^^^ = n^^\T,x\u,e^) represents the 
temperature anisotropies and polarizations caused by the metric perturbation h^, and is a function 
of the conformal time r, the comoving spatial coordinates x\ the photon frequency u, and the 
photon propagation direction e* = (sin 6^ cos 0, sin 6* sin 0, cos 6^) . 

Parallel to the Fourier expansion of hij in Eq.(l7]), n*^^) is also expanded into: 



n 
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For each Fourier component n^.^-*! 
equation is written as [4^ [54]: 



up to the first order of perturbations, the Boltzmann's 



d_ 



+ q{r)+ie'ki ) 4^^(r, u,e') 



s=l,2 



dr 



An 



(14) 



where the differential optical depth g(r) = aTNe{T)a{T) with ax being the Thomson cross 
section, and iVe(r) being the comoving number density of free electrons, /(i/) = 



dlnu ' 
s , s 

^^^^ Z]s=i '^'^'^^ hki^) being the Sachs-Wolfe term [T5] reflecting the frequency variation due 
to hij, and Vi is the velocity of scattering electrons with respect to the chosen synchronous coor- 
dinate system. In the frame associated with the density waves with a wavevector k/k = (0,0,1) 
in z direction, one has e*fj = —i^Vh with /i = cos^^ and vi, being the baryon (electron) velocity, 

. . 1 . . 2 

e^e^Pij = jj?, and e^e^Pij = {jj? — 1/3), independent of the azimuthal angle 0. Thus, as is 
expected, the Sachs-Wolfe term is ^-independent, because the k mode of density perturbation hij 
is 0-independent, as mentioned in Section 2. Then, the only term in Eq.(!T4l) that might possibly 



depend on (p is the scattering term J dQ!P{e^]e'^] 
the Chandrasekhar matrix P(e*;e'-') is only through sin( 



z'^), where the 
(h) and cosiS' - 



-dependent part of 
When one 



takes n^.^'' to be independent of 0, the 0-dependent part of J^"" d(j)'P{e^;e'^) is vanishing due to 

0. Therefore, in the case of density perturbations, it 
. This 0-independent property is the reason that the 



sm — 



/c 



cos 



Jo 

is consistent to take n^^^ independence of 
magnetic type polarization is not sourced by density perturbations (see Appendix). In contrast, 
for the case of GW, the term e^e^ Yl 



dhkjr) 



t=l,2^«J dr 



oc cos20, depending on (p. This 0-dependent 
property is also responsible for the magnetic type polarization generated by GW (see Appendix). 

To further decompose Eq. (IT4l) . one can follow the treatment of Basko and Polnarev [231 [24l [54] 
and writes n^^"* in the following form: 

/(z/)nW(z/) 



n 



(1), 



r, i^, n) 



ak{T, fi) 








(15) 



where ak is the temperature anisotropies and (3k is the polarization. By comparing Eq. (IT^ and 
Eq. (IT5l) . it is seen that, for each wavenumber k, ak is proportional to the anisotropic part of the 
intensity /, 

/(r,/i) = 7afc(r,/i) (16) 
with the factor 7 = f {v)n^'^\v) , and /3fc(/u) is related to the linear polarization Q itself 



Q(r,/i) = -ff3k{T,fi). 



(17) 
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Note that, by Thomson scattering of the unpolarized light at the last scattering, the Stokes 
parameter U = for the k mode of density perturbation. If the metric perturbation is GW, the 
form of n^.^-* will be more complicated than Eq.(IT5l). with all three Stokes parameters / = I{6,(f)), 
Q = Q(e, 0), and U = U{e, 0) ^ 0, depending on both angles 0) |221[21I2S]. See the Appendix 
for the details. 

Then Eq.([T4]) is converted into a set of two coupled first order differential equations for and 
f3k [431154], 

d \ 1 fdH ^dHi\ , , , , , , . 1 



+ q{T) + tkfij ak = -\^—- fi'^j + g(r) (^Xi(r) + tfiv, - -P2(/x)X2(r) ) , (1^ 

^ + g(r) + ikf?j h = \q{T)il - P2i^))Ur), (19) 
where P2(Ai) is the second order Legendre function, and 

1 '■^ 



2^1(7-) = - y dij,akiT,fi), (20) 

T2(r) = i j ^dfi[{l-P2{fi))^{T,fi)-P2{fi)ak{T,fi)], (21) 

play a role of sources for ak and 13^. On the right hand side of Eq.([TBj), |(^ — A*^^) is the 
Sachs-Wolfe term, which has an counterpart in the case of RGW [361 |3Zl |38]. In contrast to the 
case of RGW, the second term on the right hand side of Eq. (IT51) is a new collision term, containing 
Vb, Xi, and I2, which all contribute to the temperature anisotropies ak- From Eq.([T9l) one sees 
that X2 enters the collision term and plays the role of the source for the polarization /3k- To solve 
the set of equations ( ITSl) and ( IT9l) . one needs the quantities of g(r), H{t), Hi{t), Vbir), Ii and 
X2, which will be determined in Section 5 and Section 6. 

Although formally similar to the case of RGW [Ml EZl ES], Eqs.([TH]) and are more com- 
plicated. There are residual gauge modes contained in the solutions of H(t) and Hi(t) for the 
RD and MD eras that have to be removed before one can actually calculate and (3k- Also, the 
collision term g(r)(...) on the right hand side of Eq.f llSp is absent in the case of RGW, and needs 
some extra, proper treatments here. 

We proceed to write down the formal solutions to Eqs. (ITH|) and (11^ as the following time 
integrations 



= £rfr'e-«(-^')-^^ [1(^ - /z^^) + g(r')(Xi(r') + ^^^Vb - ^P2(^)X2(r')) 

(3k{T,fi) = 1(1 - PM) ^^rfr'g(r')e-«(-'-')-'^'=(^-^')X2(r'), (23) 



(22) 
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respectively, where 

«:(r,r')= / dT"q{T") = k{t') - k{t), (24) 



and the optical depth for the recombination 



T"0 



k{t) = k{to,t) = j dr'qir') (25) 

from the present time tq back to an earlier time r, whose time derivative yield the differential 
optical depth 

g(r) = -kir) ^ (26) 
From k{t) and g(r) follow the visibility function [401 ESI EH EZl EE ES ED] 

V{t) = -^e-"(") = g(r)e-"M. (27) 
dr 

and the exponential function 

e-^^^^ = I V{T)dT. (28) 



The quantities g(r), k{t), e^'^^'^\ and V{t) are equivalent in describing the recombination process. 
In principle, for a given cosmological model with a known a(r), once the number density of free 
electrons iVe(r) given explicitly for the detailed recombination process, one can calculate directly 
the differential optical depth g(r), k(t), and V{r) from their definitions [2Z1ES1E5]- The visibility 
function V{t) has a statistical interpretation as the probability that a CMB photon we observe was 
last scattered at an earlier time r, so that it satisfies the normalization condition 

/"TO 

/ V{T)dT = 1. (29) 
Jo 

Here we do not consider the reionization process [38], which would bring another term into the 
integrand in Eq. (!2^ . The recombination process and the corresponding V{t) depend on the matter 
fraction fi^ and the baryon fraction Q^. As depicted in FigtH V{t) is rather sharply distributed 
around the recombination time r^, and, among other things, its dependence upon the baryon 
fraction is such that a greater yields a slightly larger recombination time r^. In our context, 
by K{Td) = 1. According to [40], one has 

In practice, V{t) is often be approximated by some fitting formulae [SE HOl [371 [38] . For the purpose 
of analytic calculations of CMB polarization, V{r) was further simplified by a single Gaussian type 
of function [261 [35] 



V{r) = V{T,)exp[- \^_l^ ), (31) 
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where V{Td) = (v^vrAr^)^^ required by the normalization of Eq. (!?n) . and Ar^ is the half width 
and reflects the thickness of recombination. Ar^ also depends on the baryon fraction Qb, and a 
larger ^l^ yields a slightly narrower Ar^. It can be approximately fitted by 

Ard~ 10^^(8- 2.33 Infib). (32) 

For the redshift thickness of the recombination Az ~ 195 ± 2 by WMAPl [4j, the corresponding 
conforming time width is 2ATd ~ 0.003 fori];, = 0.045. As adopted in the case of RGW [351I371I3H]. 
for a better approximation, V{t) has also been fitted by two pieces of half Gaussian functions 

(33) 

[ \/(r,)exp(-(l^), (r>r,), 

where Ar^i = 0.0011 and Ar^z = 0.0019 for Qb = 0.045 and (Ar^i + Ard2)/2 = Ar^. It has been 
checked that the errors between Eq. (133l) and the approximate formulae proposed in Refs. [401 |5B] 
are very small, < 4% for r > r^. Eq.flH^ improves the description of visibility function by ~ 10% 
in accuracy over Eq. llHTl) . and at the same time allows an analytical calculation of C^^' . 

Now back to ak{T, fi) and f3k{T,ii) in Eqs. (!22l) and (!23l) . To get rid of their dependence of 
11 = cos9, one proceeds to expand them in terms of the Legendre functions: 

«fc(r,/i) = 5^(-z)'aKr)P/(/i), (34) 

I 

/3k{T,^i) = J2Hyi^ii^)Piif^)^ (35) 
I 

with the multipole moments given by 

21 + 1 [^^ 

aiir)=i^ — - — / diJ,ak{T,iJ,)Pi{n), (36) 



2 



1 



,2/ + l 



A(r)=^'^— / rf/i/3fc(r,/i)Pz(/i), (37) 



-1 

where the following normalization condition has been used 

^ 2 

^dxPi{x)Piix) = ^^5u,. (38) 

Inserting Eqs. (!34l) and (!35l) into Eqs. (!^ and (!?TI) . carrying out the angular integration J dfj, there, 
and using the relation (21 + l)xPi{x) = {I + l)Pi+i{x) + lPi-i{x), the sources Xi and X2 can be 
expressed in terms of the multipoles as the following: 

Xi(r) = ao(r), (39) 

X2(r)=/3o(r) + i/32(r) + ^a2(r). (40) 
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In particular, Eg. (140 ft tells that the quadrupole 0:2 of the temperature anisotropies enters X2 as a 
source for the polarization mode A(r). In the following we will express Xi(r) and X2(r) in terms 
of the scalar perturbations h^r) and i]{t). 

4. Determination of Integrands for ak{r) and Pki^) 

The Boltzmann's equation ( ITSl) and ( ITOl) can be written as hierarchical sets of equations for the 
multipole moments ai and (3i as the following. 

For each /, multiplying both sides of Eq. lITHl) by Pi{i-i) and integrating over ^ Jj*^^ d/i, one arrives 
at the hierarchical set of equations for af. 

do = -k^ai + ^h, (41) 
2 

ai = k{ao - -02) - q{ai ~Vb), (42) 
5 

2 3 2- 9 1 1 

"2 = H-ai - -as) - + 6?)) - g(r)(— «2 - -/3o - — /32), (43) 

Note that the monopole = ST/T=5y/4:, where 5^ = Sp^/p^ of the photon gas, and the dipole 
moment ai represents the velocity of photon gas. Eq. (l4TI) shows that the metric perturbation /i(r) 
induces the generation of anisotropies ao- Similar treatments of Eq. (fT9|) yield: 

/3o = -k^Pi + g(-^/3o + ^(«2 + /32)), (45) 

$i = k{f3o-'^f32)-qf3u (46) 
5 

f32 = A:(^/3i - l^s) + ?(^)(^/3o + - ^^2), (47) 

'^' = ^(^r^^'-i-4i:^^'+i)-^^'' ^^3. (48) 
As Eq.(l45il demonstrates, the quadrupole of temperature anisotropies 0^2 is the major source for 
the leading order polarization (3o via scattering. The above hierarchical sets, for both a and (3, have 
infinite number of differential equations, and should be made closed in order to find their solutions. 
One takes the cutoff 

«z = 0, A = 0, (/>3), (49) 

which is justified for the long wave modes with kr <^ 1. Dropping the small quadrupole a2, Eq. 
( 14^ reduces to 

di = kao — q{ai — ft). (50) 
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The term —q{ai—Vb) represents the momentum transfer from the baryon (electron) into the photon 
component. The quantity 1/q has the meaning of the mean free path of photons, and is a small 
parameter before the recombination. Before the recombination, photons and baryons are tightly 
coupled, and in the tight-coupling limit l/q — > 0, Eq. (15P1) implies ai = vt so that photons and 
baryons behave like a coupled single fluid [52]. But, for a more accurate account for the difference 
between photons and baryons, one keeps up to the order of l/q. To deal with v^, one needs to 
use the momentum conservation in Thomson scattering, i.e., the Euler equation for the electron 
velocity ( See Eq.(66) in Ref. [52]) 

Vb = --Vb + clPSb + ^{ai -Vb), (51) 
a K 

where R = 3pfc/4p^ is 3/4 times the baryon-photon ratio and, for a model ~ 0.045, can still be 
treated as ^ 1 during the recombination, and Cg = 1/ a/3(1 + R) is the sound speed of the photon 
gas. In the tight-coupling limit, the particle collision rate via Thomson scattering is much greater 
than the expansion rate, i.e., q ^ d/a. So, in the long wavelength limit, Eq. (!5T|) reduces to 

q 

Vb^^{ai-Vb). (52) 

Now by combination of Eqs. PTI) . (!50l) . and (!52l) . one obtains the following second order differential 
equation of the monopole: 

cio + c^k'^ao = S{t) (53) 

with the source 

S{t) ^h = H- ^-Hu (54) 

where the second equal sign follows by Eqs. lITTl) . Note that appearing in Eq. (l53l) is a function of 
time through the ratio R. At the level of our analytical calculation, R will be approximately treated 
as a constant, and its value will be taken at 5; ~ 1100 around the recombination. The general 
solution of Eq.(!53|) is given by the following simple form 



ao(r) = B, cos(j9r) + sin(pr) + f '^(^0 sm(pr pr') ^ 

Jo V 



(55) 



where p = Cgk, and Bi and B2 are constants, to be fixed by initial conditions. As we shall see 
later in Section 6, both Bi and B2 can be set to vanish. Once the scalar perturbation mode /i(r) 
is specified, the monopole «o(t) follows from Eq. (l55l) . and so does the source Xi(r) via Eq. (!391) . 
We remark that if the expansion term —^Vb in Eq. (l5ip was kept, Eq.flSSD would be modified as the 
following 

R R \ ■ 

"0 + Y^-^tto + c^/c^ao = Y^-^-/i + 5(r). (56) 
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With R being a time-dependent function, this differential equation could also be solved, whose 
solution would differ only slightly from Eg. (1551) on large scales under consideration (/ < 400). In 
order to get full-analytical formulae, we use Eq. (!55l) . 

Once ao^r) and /i(r) are known, the dipole ai follows immediately from Eq.(|4l 

/ N 3, . 1 

Let us calculate the source X2 for polarization in Eq. (l4DI) . It is interesting that a linear combination 
of Eqs. (l43l) . (l45l) . and ( l47l) . with higher order terms (/ > 3) in perturbations being dropped, leads 
to the following differential equation 



-h). 



(57) 



i2 + = M(r) 



(5^ 



with 



M(r) = --Mr) - ^V{r) = -^«o + ^H. (59) 

Eqs.( l58l) and (14U1) . tell that M(r) is the source of (3q and of a2, simultaneously, and is expected 
to contribute equally to them as well. Eq.(!58l) has a formal solution 



X2(r) = r dr' M(r')e-^"("'"') 
Jo 



(60) 



with k(t,t') being defined in Eq.(!24l). As will be seen later, M(r) is basically contributed by 
the gradient of the peculiar velocity of photon fluid, kai^r). When the perturbation mode r]{r) 
is specified, one calculates I2 from Eq. (!5Dl) straightforwardly. Having obtained Xi and I2, one 
proceeds to perform the time integrations of the modes and (3^. 

5. Time Integrals for Temperature Anisotropies and Polarization 



As demonstrated in Appendix, by projecting and /3k on the basis Pi{fi) in Eqs. (ll40p and 
(1141 p . one obtains the multipole moments of and af for the temperature anisotropies and the 
electric type of polarization respectively, which have been given by Eqs. (ll33p and (I139p as the 
following 



dr 



d^ 



d 



d^ 



+ i/.^Je-^^^) + n.)^X, + -X.(l + - 



3 /(/ + 2)! 



1/2 



dr V{t)X2{t 



MO 



4 V(/-2)!; Jo V ^ - ^ ^2 ' 

where the variable ( = A;(ro — r). These two time integrations have to be carried out. 
First, we calculate af. Substituting I2 of Eq.(l60l) into Eq.(!62l) gives 

3 / (/ + 2)! ^ f-^" 



(61) 
(62) 



drV{T)^^ r dT'M{T')e~^o-(-')+w-(-). 
C Jo 



(63) 
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As we have seen, the visibility function V{t), as an integrand, is narrowly peaked around the 
recombination time with a width Ar^, so, for the r-integration in Eg. (1631). the integrand around 
Td only will have significant contributions. Moreover, in the r'-integration the exponential factor 

Q-jo'^i'^') behaves like a step function: e"^''^"^'-' ~ for r < r^, and e"^'^'-^ '' ~ 1 for r > r^. Thus, 
the integrand factor M(t') can be approximately pulled out of the r'-integration, leading to 

- i (f^) f ''^"(^'f^^' W [ e-**'.-**>..'. (64) 

To perform the r'-integration in Eg. (1640 . one introduces a variable x = -^^^ to replace the variable 
r' im [2H]. The corresponding limits of integration are r' = r — ^ x = 1 and r' = — )■ x = oo. 
Since V{t) is peaked around r^ with a width Ar^, one can take dr' — ^Ar^ as an approximation, 
valid over the period Ar^ around the recombination. For a justification of this approximation, in 
Figl2]we plot the optical depth k(t) as given in RefflD], which indeed behaves approximately as 
an exponential: k(t) oc e^^/^'^'* around r^. Then 

. I (|±|) " A., £ *^,.,f)M(., e-*.--*^W. (05, 

V{t) contains e~'^^'^~'^'^^' , and j;(C) contains a mixture of oscillating modes e*^"^ and e"*^^ with 
p (X k. Using the formula of a form e~'^'^^e*^^(ir = e^^^ e~'''^^(ir, the r-integration is 
rendered approximately into 

r rfry(r)^M(r) ^ M|^M(r,)D^(A;) T rfrV(r), (66) 

where 

is the Silk damping factor [62] for the polarization. It arises because the CMB photons diffuse 
through baryons and smaller scale fluctuations are smoothed, i.e., those modes of higher k are 
more effectively suppressed. Mathematically, it occurs as a sort of the Fourier transformation of 
V{t) in Eg.(l33l). This is one of the advantages of our calculation in that the Silk damping factor 
arises naturally instead of adding by hand, ce and Be are two fitting parameters, and the values 
ce ~ 0.27 and Be ~ 2.0 yield an agreeing match with the numerical results [111 [T5] over a range 
I < 500. The physical interpretation of the appearance of D(k) associated with the recombination 
process has been given in Refs. [351 EZ]- During the recombination around the time r^, the last 
scattering of CMB photons off baryons occur effectively only within a time interval ~ Ar^. Putting 
it in terms of the spatial scale, the smoothing of density fluctuations by the associated diffusion 
through baryons occur effectively only on a scale of the thickness of the last scattering surface 
~ Arrf (note that we use unit c = 1). Those modes, e*''"^ and e^**^^, with wavelengths shorter than 



15 



~ At^ are effectively damped, whereas the long-wavelength modes are less damped. We remark 
that, as a fitting formula, DE{k) in Eq.(!67l) works only approximately, since other time-dependent 
factors in the integrand have been treated as constants. Besides, there are other processes [121122], 
which are significant on small scales, are not taken into account here. So the parameters ce and 
Be are introduced in Eg. (1670 for adjustments. Among the two terms in Eq. (l67p . DE{k) is more 
sensitive to the term with a smaller time interval Ari. 

The remaining double integration in Eq. (!65|) can be carried straightforwardly 

r dTV{r) r — e-^^i'^W^+i^'^W = r dK e-^^^ H — e"^^^ = ^ In — , (68) 
Jo Ji ^ Jo Ji ^ 7 3 

whereby Eq. (!?n) has been used in the first equality. The above treatment of the integrations is 
similar to that in Refs.[3ni [3Zl ES] for the case of RGW as the source. One arrives at the explicit, 
analytical formula of the multipole moment of polarization 

which depends upon the function Mij^) at the recombination time t^. As a marked feature, af 
contains explicitly the recombination width Ar^, which arises from the r'-integration in Eq. (!54f) . 
Since Ar^ is small, the amplitude of af will be consequently small, in comparison with aj , whose 
dominant part does not contain this Ar^ as will be seen later in this Section. Physically, the 
factor M(r) represents the source of both the leading order polarization (3^ and the quadrupole 
temperature anisotropies a2, and contributes equally to them as well. As its time accumulated 
effect, the factor ATdM^Td) appears in af in Eq. (!5^ and in the last term of af in Eq.lJTHI) as 
the quadrupole temperature anisotropies. During the course of time, the contribution of M(t) is 
significant only around the recombination time Td with a width Ard- Note that the spherical Bessel 
functions ji{k{To — Td)) in Eq. (l69p is narrowly peaked around A;(ro — Td) — I for / ^ 1. In our 
notation tq — ~ 0.97. So, for each given multipole /, the factor ji{k{To — Td)) serves as a filter, 
selecting those modes with a wavenumber A; ~ / for af. 

Next, we calculate af in Eq. (!Ml) . The first term in the integrand of Eq. (l6ip is the integrated 
Sachs-Wolfe (ISW) contribution, and contains the exponential factor e~'^^'^\ which can be roughly 
approximated by a step function [371 135]: 

e-W = I °' ^< J (70) 

1, for Trf < r < To. ^ ^ 

So the ISW term is approximated by 

£ dr (h + H,-^^ e--MjKC) = dT (hMQ + i^.^jKo) , (71) 
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where the lower limit has also been replaced by r^. In the pertinent domain, H{t) and Hi{t) are 
comparable to each in magnitude, whereas in Eq.(!TlIl. the integrated value of ^j/(C) is two orders 
of magnitude smaller than that of so the term ^j/(C) in Eq. (!7ll) can be neglected in the 

estimation. So the left-hand side of Eq. (!7ll) reduces to 



drHjiiC) - (Hiro) - H{Td))3i{kT,). (72) 
The second integration in Eq. (!6TI) is 

jyrV(r) L + a.|-f £e-*"<-'.AW (l + 



m, (73) 



where we have substituted Xi = a^, vi, = ai in the tight-coupling limit, and X2 as given in Eq 
The first two terms in Eq. (!73l) can be treated as before, yielding 

dTV{r)ao{r)ji{0 - a,{T,)DT{k)ji{k{T, - r,)), (74) 

r° d d 

/ dTViT)a^{T)-JliC) ^ a,iT,)DTik)-UkiT, - r,)), (75) 

Jo «C "C 

where the damping factor for the temperature anisotropies is taken as 

with ct and bx being two fitting parameters, and ~ 0.65 and bx ~ 0.6 yield a good match 
with numerical results by CAMB over a range / < 500. The last term in Eq. (l73p is a double time 
integration and has the same structure as af in Eq.( 15^ . and can be treated in the same way, 
yielding 



Q /"TO / J2 

dTV{T)(l + ^]jl{C) I e-^o-(ry)M{T')dT' 



-l^ln^Ar,M(r,)D^(fc) (l + ^) j,(C)lc= 



^fc(fO-Td)- 



(77) 



Again, the term proportional to -^ji{C) can be neglect when calculating the power spectrum in 
Section 7. Putting these four pieces together, one arrives at the explicit, analytical formula of the 
multipole moment of temperature anisotropies 

aj = ao{Td)DT{k)ii{k{TQ - Td)) + ai{Td)DTik)^ji{k{To - Td) 

+ {H{to) - HiTdmkro) 

15 10 f d"^ \ 

--ln-Ar,M(r,)D^(fc) [l + — j 3m\c=Kr.-r,)- (78) 

In the above expression, the term is dominant, the ai term is secondary, the ISW term is smaller 
than the ai term, and the last term containing the factor Ar^^ is smaller than the ISW term. 
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Note that the two major terms, oq and ai, in Eq. (!78l) do not contain Ard- This is because their 
corresponding integrations, Eqs.(!74D and ( iTSf ). are single time integrations, instead of double time 
integration. By comparison, the amplitude of af is expected to be higher than that of af . We 
note that the structure of af in Eg. (1751) is similar to the parallel formula in the Newtonian gauge 
given in Ref.[40j, which did not have the last term oc ATdM^Td). The relative contributions of the 
four terms in Eg. (1751) will be demonstrated in FigEl 

To completely determine af in Eg. (!69l) and af in Eg. (!78l) . one still needs ao, ai, M(r) in 
Egs.(l55l). ( 1571) . and (15^ . respectively, which all depend upon the scalar perturbations /i(r) and 
?7(r). In the following we will solve for /i(r) and 77(t). 

6. Determination of scalar perturbations 

The unperturbed spacetime background are described by the Friedmann eguations: 

-Ga'^p, (79) 



a 3 
fi 4-7r — 

-= Ga'ip + 3P), (80) 
a 6 

where p and P are the mean energy density and pressure. The Einstein eguations for the scalar 
perturbations in synchronous gauge are the following [491 [52] 

er]-^^h = 47rGa^6T^, (81) 

k^T] = AnGa'^iySTf, (82) 

h + 2-h- 2^7] = -MGa^5Tl (83) 
a 

h + 6ii + 2- (h + 6f]^ - 2^7] = -2ATiGa^{p + P)(t, (84) 

where o represents the anisotropic stress, 5Tq = —p5 is the perturbed energy density, 5T- = 35P 
is the perturbed pressure, and ^ = c^, where the sound speed — | in RD era and — in MD 
era. 

First, let us do for the RD era. We are concerned with the long wave modes with kr <^ 1, and 
the solutions of the set of Egs.(l8l])-(|84D for /i and are [47l l48l l49l l52] 

h = A + B{kT)-^ + C{kT)^ + D{kT), (85) 
7] = 2G + ^D{kT)-\ (86) 

All the coefficients A through D actually depend on the comoving wavenumber k, which has been 
skipped hereafter for notational simplicity. The two terms proportional to A and B are gauge 
modes, which will be dropped, and two physical modes are proportional to C and D. 

h = G{kTf + D{kT). (87) 
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Among these two modes, the mode (fcr)^ grows faster and is dominant at late times, and the 
mode (kr) is less important, which was neglected in the treatment of Ref.[52j and was taken to 
be small in Ref.[49j. In principle, the two coefficients C and D should be determined by either the 
inflationary or the reheating era that precedes the RD era. To avoid further complication from the 
preceding eras, we shall treat D as a small parameter proportional to C. For simplicity of analytical 
calculations, we do not include the modifications due to cosmic neutrinos, which will bring higher 
order terms (kr)'^ to r] during the RD era |52]. Let us examine the long wave approximation during 
the RD era. At the radiation-matter equality r2 the comoving sound horizon is ~ CsT2. Those 
/c-modes with 1/k > CsTi can be taken as the long wave modes during the RD era. In our notation 
with the comoving time r specified from Eq.([T]) through Eq.(l3]), this is equivalent to A; < 210 (0.025 
Mpc~^). For wave number greater than this, a more elaborated treatment of the perturbations 
during the RD would be desired than presented here. 

For the MD era, the solution of the metric perturbations are given by |471 |4H1 US |52] 

, = + _L_f+_L_G, (88) 

„ = 5£ - j^,F, (89) 

where the constant J is a gauge mode corresponding to a transformation of the spatial coordinates, 
i.e., a rescaling of the scale factor a(r), and can be dropped. As has been known [47], for /i(r) in 
Eq. (!88l) . the linear combination + -j^^ is another gauge mode, which is dominated by -f^p^ for 
A;r ^ 1 (long wavelengths or early time), and by for kr \ (short wavelengths or late time). 
In our context, we aim at the large angular temperature anisotropies and polarization of CMB. So 
we are concerned with the long wavelength perturbations around the radiation-matter equality T2 
and the recombination time r^. Thus the G term is taken as the dominant gauge mode. To keep 
our analytical calculation simple, we drop the G term. In fact, the G term is the time-translation- 
invariant solution and can be gauged away by a restricted coordinate transformation within in the 
synchronous gauge [4H]. Other discussions on gauge modes are given in Refs. [471 [45]. The term 
proportional to E in Eq. (!88l) grows with time and is the primary portion of the physical mode. 
Thus, for the MD era, one has 

h = (krfE + -i-F, (90) 

'/ = 8£ - (^f- (91) 

From these specifications, the source S{t) in Eq.(l54l) reduces to 

2A-2 

S{t) = C, for RD, (92) 

3 
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91.2 2 

S{t) = —E + ^F, for MD. (93) 

Now we need to make a proper connection of the perturbations for the RD and MD eras at the 
radiation-matter equality r = T2. We remark that the energy density perturbation S is continuous 
in the transition from RD to MD era. But the pressure P is not required to be so, as P > 
during RD, and P = during MD. By the perturbed Einstein equation Eq. (!5T]) for 5, one finds 
that the combination k'^rj — is required to be continuous at r = T2. Since d(r) is continuous 
as prescribed in Eqs.(I2]) and ([3]), h and rj are required to be continuous, leading to 



p 

2k'^T2C + kD = 2A;VE ^ , (94) 

kr^ 



From these two algebraic equations, one solves for E and F in terms of and C and D. 

F=-(kT,f(2C+J^) . (97) 



QkT2 J 

The coefficient D is a small parameter that needs to be fixed. We take the coefficient D to be 
smaller than C by a factor (A;r2) in the long wavelength limit kr <^ 1, so that D ~ {k,T2)C . 
Specifically, in the following analytical calculations, we take 

D = -y (A;r2)C, (98) 

though other possible choices may also be justified as long as D is subdominant to C in the long 
wavelength limit. Substituting Eq. (!9Hf) into the above yields 

E = \c. F=^^{kT2fC. (99) 

One can check that, in the RD, as well as in MD era, if we transform the perturbations h and rj in 
synchronous gauge back to the and if) in Newtonian gauge [4D], the results are consistent with 
each other. Fig. |4] shows the continuous joining of the perturbation modes hij) and rfij) at T2, 
and Fig. [5] shows the continuous joining of the modes H{t) and Hi{t). As one can check, the 
functions h and fj are not continuous at r = r2 by our joining condition. 

To fix the initial condition, we need to specify the fc-dependent coefficient C. According to 
the inflationary models of the early universe, the primordial scalar perturbations were generated 
with a nearly scale-invariant spectrum with a spectral index ~ 1 [64j- In our notation this 
corresponds to C oc A;^'^"^"^^ For inflationary models proposed so far, the most uncertain quantity 
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is the amplitude of the spectrum. In practice, this can be fixed by cosmological observations, say, 
the WMAP result. One writes the curvature perturbation spectrum 

/ 7 \ ns-l+|Qs ln(fc//co) 



where the physical pivot wavenumber fco = 0.002Mpc ^, A^(A;o) is the normalization at /cq. 
WMAP5 [7J gives A%{ko) = (2.41 ±0.11) x 10"^ WMAP5+BA0+SN Mean [9J gives Al{ko) = 
(2.445 ± 0.096) x 10^^. Besides the scalar spectral index Ug, we include a possible scalar running 
spectral index in the spectrum |1T1I22]- The fitted value of is much affected by the presence 
of as and the RGW component, and by additional combination with SN la and BAO data as well. 
In absence of a, and the RGW, WMAP5 gives = 0.9631^1^^^ [7], WMAP5 +SN la+BAO gives 
n, = 0.960l|J:S}^ and WMAP7 gives = 0.963±0.012 [9J. When a, is allowed, WMAP5 gives 
Ug = l.OSltooll and as = — 0.050±0.034 with a better determination of the third acoustic peak 
|9], WMAP5+BA0+SN has given = 1.089;|J;|J^^ and as = -0.053+|]:[J^^ |Zl ^. More recent 
WMAP7+ACBAR+QUaD gives = 1.041^°°^^ and a, = -0.041+°:°| When RGW is also 

allowed p4|, WMAP7+Tensor gives ra, = 1.076 ±0.065, a, = -0.048 ± 0.029, r < 0.49 [9j. In the 
slow-roll scalar inflationary models, Us and can be calculated from the inflationary potential and 
its derivatives [1T1[12J. For generality, we will treat and as parameters. Eq. fllOOp corresponds 



where the normalization Co ~ 204. The physical pivot wavenumber ko corresponds to a comoving 
wavenumber kc = fcoa(To) — 17.1 for the Hubble parameter Hq = 70.1 km s^^ Mpc~^ [9]. 

In fixing the initial condition for ao at kr <^ 1 during the RD era, the coefficients Bi and 
B2 in Eq.(l55l) have to be specified. In the tight-coupling approximation, kai in Eqs. (l4T]) can be 
neglected, yielding ao = h/3 for kr ^ 1. By comparison, in the limit kr ^ 0, h behaves as in 
Eq. (157l) . so the term Bicos^pr) in should be vanishing, leading to Bi = 0. The B2sm{pr) 
term in Eq. (!55l) represents the isocurvature mode of initial perturbations. A stringent constraint 
has been given by WMAP5 on the isocurvature contribution with the isocurvature/adiabatic ratio 
a-i < 0.015 at 95% CL |71[38]- For simplicity, we can choose the coefficient B2 = 0. Then the 
monopole ao{r) in Eq.(!55l) reduces to the integration 




(100) 



to 




(101) 




dr', 



(102) 



Using Eqs.(IS2I), (ES]) into the above yields the monopole 

[3 2 
aoi^r) = 2C(1 + R) - cos(pr - PT2) - cospr + - 




(103) 
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From Eq. (!57l) and Eq 



follows the dipole 



«i(r) =2C^'i{l + R) 



- sin(pr — PT2) — sin(pr) 
5 



.36 



(104) 



By the definition in Eq. (l59p . we take time derivatives of ao{T) and ?7(r) in Eq. (11031) and Eq. 
respectively, and arrive at 



M(r) 



--Cil + R) 



— p sm( pr — PT2) + p sin pr 
5 



72 
— ( 

25 



C(l + i?) / ^ cos(pr - pr')d(pr') - C 

J TO 



216 r| 
"25" 74 



(105) 



for the MD era. We have checked that, in this final expression, M(r) is dominated by the first 

3 

two terms coming from — 2q;o(t)/5, whereas the last term —C^^ coming from — 4?7(r)/5 is 
comparatively small by more one order of magnitude. It is important to notice that, due to time 
differentiation, M(r) contains the functions like sin(pr—pr2), whereas aQ^r) oc cos {pr —pr2) + .... 
This fact will lead to the character of the present CMB that the peaks of the polarization Cf"^ oc 



,E|2 



OC I sm(pTci — ]5T2)p + ... and of the temperature anisotropies Cj^-^ oc |af P oc | cos{pTd 



PT2)P + ••• appear alternatingly. 

Using Eqs.dm, (jSU]), (EI]), and 
MD, are given by 



, the time derivatives of the scalar modes, H and Hi during 

1 0S 

H{r) = (106) 



H,(t) 



(107) 



('2A;V - 9^4 + 1624 
' 5 \ r"' 

So H{r) and Hi{r) are comparable to each in magnitude around the recombination. 

7. The analytical spectra 



Given the multipole moments af in Eq. (l69|) and af in Eq. (!78l) . the spectra Cf^, Cf^, and 
Cf^ are calculated as the following integrations over the wavenumber k 



TT 



\a[{k)\^kdk, 

a'[(k)af{k)kdk, 

\af{k)\^kdk. 



(108) 
(109) 
(110) 
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The resulting spectra are explained in the following graphs. 

Fig. |3] demonstrates the relative contributions by each term to Cf^ . Over the relevant range 
I < 500, the contribution by ISW is rather flat as a function of /, and its amplitude is at most 
~ 10% that of the ai term. The last term in Eq. (!7Hl) contains the factor /^TdMir^) given by 
setting T = Td \n Eq. dlOSp . and its contribution to aj is even smaller than the ISW term, with two 
low bumps at / ~ 130 and at / ~ 350. The smallness of this term is due to the extra small factor 
Ar^ ~ 0.003. Thus, the major features of aJ in Eq. (!78l) are largely contributed by a;o(Td) and 
Cii{Td), whereas the quadrupole part of aJ contains the factor ArdM^Td), similar to the polarization 
af in Eq. (16911. This also tells that the polarization af is smaller than the temperature anisotropies 
af in amplitude. 

In FigHwe plot these analytical spectra /(/ + 1)C[^ , 1(1 + l)Cf'^, l{l + l)Cf^ for n, = 0.96, 
tts = 0, and the baryon fraction ^l^, = 0.045. For comparison, the numerical result from CAMB 
[T5] and the observed result from WMAP5 [7] are also given. For a more realistic case, one would 
have to also include the analytical C^^' by RGW [361 |3Zl |38J at a tensor/scalar ratio r to form the 
complete calculated Cf^'. We leave that for future studies. Fig. [7] shows that the overall profiles 
of the analytical spectra agree well with the numerical and the observed on large angular scales with 
/ < 500. This range covers the first primary peak of Cf"^ and the first two primary peaks of Cf^ 
and of Cf'^. Only around / ~ 310 where the second primary peak of Cf^ is located, the analytical 
Cf'^ deviates by ~ 18% higher in amplitude from the numerical one. For smaller angular scales, 
the analytical results deviate considerably from the numerical ones. This has been expected since 
our calculation is based upon the long wavelength approximation valid only for large angular scales. 
From Figj6l we see that the first two peaks of Cf^-^ occur at / ~ 200 and / ~ 500, while those of 
Cf^ occur at / ~ 100 and / ~ 400. This alternating occurrence of the peak locations of Cf^ and 
Cf^ has been anticipated. (See the discussion below Eq. dlOSp . Based on the analytic results, one 
can estimate the span of the two adjacent peaks of Cf^ in /— space, which corresponds to that 
of I sm(csk{Td — T2))p in k— space. Since ji{k{TQ — Td)) is significantly contributive only around 
A;(ro — Trf) ~ / for / ^ 1, it plays a role of a filter and selects those A;(ro — Td) ^ I part of the 
integrand to contribute to the integration J dk over k. Qualitatively, the span AA; of two adjacent 
peaks of | sm.{csk{Td — T2))P is given by a relation tt = CgAk^Td — T2). Then the span of the two 
adjacent peaks of Cf^ in /— space is Al ~ Ak{To — Td) ~ 370. The same Al holds also for Cf^. 
This is roughly what is seen in FigH (See also Ref.[65]). 

In FiglZl we sketch the profile of /(/ + l)Cf^^ as a function of/, which, notably, has two bumps, 
one at / ~ 100, and another at / ~ 400. In order to interpret the origin of these two bumps, we 
also sketch the main factor DEik)M{Td) /[k(To — Td)]"^ of af in Eq.(!69l) as a function of A;. By the 
projection of j/(A;(ro — Td)), the square of DEik)M{Td) /[k{TQ — Td)]"^ around k I, aside some 
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factor, is basically Cf^ around /. Since DE{k)M{Tci) /[k{TQ — Tii)]'^ has two bumps, around k ~ 100 
and k ~ 400, they give rise to the two bumps of Cf"^. 

FiglH] shows the first two peaks of the squared time derivative fc|do(T)p. Below Eq. fllOSp we 
have mentioned that — 2do(T)/5 is the dominant term of M(r). Thus, it is the time derivative 
aoi^r) that essentially determines the characteristic profile of C^^^, including the peak locations. 
The two peaks of |do(T)p consequently gives rise to the first two peaks of C^^, the first one 
actually being very low so that it is only a low bump. 

f \gM shows the dependence of Cf^' upon the scalar spectral index n^. The pivot point 
fco = 0.002Mpc^^ corresponds to / ~ 12. As is seen, a greater value of yields a higher amplitude 
of C^^' for / > 12. This is expected from the initial amplitude C given in Eq. dlOip . which gets 
larger for a greater Ug in the range k > k^c- The effect is most obvious around the primary peaks. 

FigJTU] shows the dependence of Cf-^' upon the scalar running spectral index a^. A greater 
yields a higher amplitude of C^-^' as is expected from Eq. fllOip in the range k > k^c- Comparing 
Figj9]with FigJTO] reveals that there is a certain degree of degeneracy between the indices Ug and 
Us as two major cosmological parameters. This degeneracy has demonstrated itself in fitting the 
observational data of WMAP [6l |71 El |9] . Therefore, given the accuracy of current observational 
data of C^^' , it is not easy to distinguish the fine details of the inflation potentials. 

FigHU shows the dependence of C^^' upon the baryon fraction ^Ib, in the amplitudes and 
the locations of peaks and troughs. As is seen, a greater value of ^l^ yields higher amplitudes 
of Cl^ (also see Refs. [4U],|55], [55], [5B]) and C[^, but a lower amplitude of Cf^. This can 
be understood as follows. Eqs. (ll03p and (11040 show that a greater corresponds to a greater 
R and gives higher amplitudes of ao and ai, hence a higher amplitude of af in Eq. (!78l) and of 
Cf^. On the other hand, af in Eq.(l69l) is proportional to the recombination width Ar^, which 
is smaller for a greater value of Vlh as fitted by Eq.(!32l). Thus Cf^ has a lower amplitude for a 
greater As for the cross spectrum Cf^^, the fi^-dependence of its amplitude is the outcome 
of these two competing factors. Since greater Ug and Vlb both tend to enhance the amplitudes of 
the spectra Cf^^ and Cf^^ , there is also a degeneracy between and Vlb in regards to Cf^ and 
Cf'^. Nevertheless, for the spectrum Cf^, greater Ug and have just opposite effects on its 
amplitude. This feature will help to break the degeneracy. Fig{TT]also shows that a greater f2b 
shifts the locations of peaks and troughs of C^^' to larger / (smaller angles). This is because a 
greater fi^ leads to a lower sound speed Cg of photon gas, so at a fixed frequency the corresponding 
wavelength is suppressed |66j|. By the analytic results, this is evident from the oscillating factors 
sin(csA;ro) and cos(csA;ro) contained in af and af , whose peak locations are stretched to a larger 
wavenumber k ( i.e., larger / via the projection of j^(A;(ro — r^)) ) for a smaller Cg. 

FigJTO shows that a longer recombination process (a greater Ar^) yields a higher amplitude of 
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polarization. This property has been obvious since the analytic expression in Eg. (15^ tells af oc Ar^. 
FiglT2]also shows that a longer recombination process brings more damping of Cf^ on small scales. 
This is because af in Eg. (15^ contains the damping factor DE{k) oc e~'^^*^'^^'^'*^^ Similarly, this 
feature also is shared by Cf^, as the damping factor DE{k) oc e~'^'^^^^'^''^^ appears in the major 
term of af in Eg. (El]). 

FigJT^ shows that a longer Ar^ yields higher peaks as well as lower troughs of Cf^. Moreover, 
a longer Ar^ slightly shifts the peaks and troughs to larger scales and causes more damping on 
smaller scales. These features are helpful to probe Ar^, as long as current and future CMB obser- 
vational data are accurate enough. However, as an approximation, this analytic result also has its 
limitation, since the recombination history has been primarily represented by only two parameters: 
the recombination time and recombination width Ara as an integrated effect. Two different re- 
combination histories via different differential optical depth g(r) would lead to the same amplitudes 
of bumps and troughs, as long as they have same and Ar^^. 

FigjM] shows that a late recombination (greater r^) shifts the peaks and troughs of the po- 
larization Cf^ to larger angular scales. The property also holds for Cf^ and Cf^. This can be 
explained by the appearance of the function ji{k{TQ — r^)) as a factor in the analytic expressions of 
af and af . 

FigJTSl shows the ratios of the analytic spectra to the numerical spectra, Cf'^{a)/Cf^{n), and 
Cf^{a)/Cf^{n). The ratios are seen to be centered around 1 for / < 500, showing a reasonable 
agreement between the analytic and numeric on large angular scales. 

8. Conclusion and Discussions 

In this paper, we have presented an analytical calculation of CMB anisotropies and polarization 
generated by scalar metric perturbations in the synchronous gauge, resulting in the explicit, analytic 
expressions of the multipole moments af in Eg. (178 p and af in Eg. (15^. and, thus, of all the analytical 
spectra Cf^, Cf^, and C™. This has been implemented primarily through an approximation 
treatment of time-integrations over the recombination process, a technigue used before for the 
case with RGW as the generating source [361 [371 ES]. We have also dealt with the removal of the 
residual gauge modes and the joining condition at the eguality of radiation-matter of the scalar 
perturbations. Several approximations have been used, such as the long wavelength approximation 
for scalar perturbations during the RD era, tight-coupling approximation for the photons during the 
recombination process. 

These results are new and have significantly extended the earlier preliminary works. The analytic 
expressions of polarization af and the related spectra Cf^, and Cf^ are what have not been 
addressed in Ref.[42j. Besides, our analytic expression af fulfils what was not completed in Ref . [43] . 
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and, to a great extent, improves what was given in Ref.[45], as our expression af contains the 
separate contributions of monopole, dipole, quadrupole, and Sachs-Wolfe terms. 

Our analytic calculation shows that the polarization af is generated mainly by the quadrupole 
of temperature anisotropies a2 via scattering. Besides, af and of ^2 are simultaneously generated 
by the combination M(rd), so that the resulting af and 0:2 have a similar structure and both are 
smaller than the total temperature anisotropies af . 

Furthermore, the analytic expressions of af and af demonstrate explicitly that the peaks of the 
polarization Cf^ and of the temperature anisotropies Cf^ in /— space appear alternatingly. These 
help to understand the important features Cf^' . 

As the major advantage of analytic expressions, af and af explicitly show the dependance upon 
the scalar perturbations, initial amplitude Cq, primordial spectrum index Ug, baryon fraction ^Ib, 
damping factor D{k), recombination width Ara, and the recombination time Td- These properties 
are transparent in analytic expressions, but might not be directly obvious in the numerical code 
itself. For instance, the dependencies upon Ar^ tell that a longer recombination process yields a 
higher amplitude of polarization since af oc Ar^, and brings more damping of af and af on small 
scales through DE{k), D'r{k). The dependencies upon tell that a late recombination shifts the 
peaks and troughs of spectra Cf--^' to larger angular scales. 

The spectra Cf^, Cf^ and Cf^ agree with the results of the numerical codes on large angular 
scales / < 500, covering the first two peaks and troughs of Cf^'. On smaller scales, the analyt- 
ical spectra deviate considerably from the numerical ones, as is expected for the long wavelength 
approximations. Serving as a complement to the numerical studies, the preliminary analytical cal- 
culations efficiently promote the analysis of effects upon Cf^' by various physical processes, and 
improve our understanding the important features of the observed CMB. 

Based upon the framework presented in this paper, several points can be further improved for 
more accurate spectra Cf^'. Some of them are listed as the following. One can extend the 
analytical calculation to smaller scales [SZ]- For the solutions of perturbations h and r], Eqs. (!55l) 
and (1551) . one may include higher order terms in kr. Consistent with this, one could do a finer 
treatment of the baryon component before the recombination, including the time dependence of 
the ratio R{t) as in Eq. (!5^ . One could also try to include the modifications from the relativistic 
neutrino component during the RD era. Finer examinations can be made on the initial condition 
during the RD era. For instance, alternative forms could be tried for the slowly growing mode 
D other than that in Eq. (!^ . and possible allowances could be tested for initial isocurvature 
perturbations besides the adiabatic ones. Further examinations on the gauge modes for smaller 
scales could be made during the MD era. Very importantly, one should include the reionization 
occurred around a redshift z ~ 11, a process secondary only to the recombination. This will 



26 



definitely bring about modifications of C^-^' on large angular scales / ~ 5 [3B]. Finally, to extract 
possible signals of RGW from observations, one should separate the contribution of RGW with 
various ratio r from scalar perturbations in the total spectra Cf^', which can be done within the 
framework in synchronous gauge by using the results in this paper and our previous work on RGW. 



Appendix: The multipole moments for radiation field 

On a 2-dimensional unit sphere with a metric 

da^ = g^f^dx^dx^ = de^ + sin^ ed(f)\ (111) 

a general radiation field is usually characterized by the following 2x2 polarization tensor [531 l29l l43] . 

p _if I + Q -{U-tV)sm9\ 

with the four Stokes parameters (/, Q, U, V), where / is the intensity of radiation, Q and 
U describe the linear polarization, and V is the circular polarization. In the case of CMB, the 
Thomson scattering during the recombination does not generate the circular polarization [53J, so 
we set V = 0. Note that / is a scalar on the 2-dim sphere under the transformation of 9 and 0, but 
Q and U transform among themselves. To deal with this problem, several formulations have been 
proposed, such as the total angular momentum method using the spin-weighted spherical harmonic 
functions |68], and the spin raising and lowering operator method [281l29j. These two treatments 
are essentially equivalent, and the latter will be adopted in the following. The tensor in Eq. (lll2p 
consists of two parts: 

Pab = -iQab + Pab^i 

where \lgab foi' the temperature anisotropies is of scalar nature, and -Pj,^^ for the polarization is 
the symmetric trace-free (STF): 

^\-Usm9 -Q 8111^9 J' ^ ' 

from which one can construct two linear independent, invariant quantities involving its second order 
covariant derivatives [43J: 

Ei9, <P) = -2{P!ry^^ B{9, <P) = -2{P!ry'-\l (114) 

where 

/ -sin6' \ , - 
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is a completely antisymmetric pseudo-tensor. £" is a scalar on the 2-sphere and i? is a pseudo-scalar. 
It is revealing to write Pa = (Pab^^)'^- Then E = Pa '"^ '\s a divergence of Pa, and B = Pa-b^"'^ is a 
curl of Pa- In this regard, E is referred to as the "electric" polarization, and B as the "magnetic" 
polarization. As can be checked, by directly calculating the covariant derivatives on the 2-sphere, 
one has _^ 

E = -\[^ \Q + iU)+d \Q - lU)], (116) 



B = -[d\Q + iU)-d\Q~iU)l 
where d'^ \s the raising operator acting twice, and 8^ \s the lowering operator acting twice. 



{Q-iU){fi,(f)) = {-d, 



1 — /i' 
-id. 



:)1(l-/x^)(Q-^t/)(/i,, 



d'{Q + ^U){^^, 0) = {-d, + ^^mi - /x^)(g + ^t/)(/x 



(117) 

(118) 
(119) 



where fi = cos 6. 

Since /, E, and B are scalars on the 2-sphere, they can be expanded in terms of the spherical 
harmonics Yim{0,(l)) as a complete and orthonormal basis 



(120) 




(121) 



(122) 

I /, — Z I ! 

1=2 m=-l 

For technical simplicity, one can choose the coordinate with the polar axis z pointing along the 
wave vector k of the scalar perturbation mode: k||z. Let an unpolarized incident light have 
an intensity /' scattered on a charge. Using the differential Thomson scattering cross-section 



^ = ■ e|. with e' and e being the polarization of incident and outgoing light, respectively. 



3i7t I 



Stt 



one obtains [27j / = ^{6)1' {I + cos^ d), Q = ^{9)1' sin^ 9, and f/ = for the outgoing wave, 
where 9 is the angle between the incident and outgoing directions. The result does not depend on 
the azimuthal angle (p. As a corollary, in an azimuthal symmetric configuration, Thomson scattering 
of an unpolarized light yields 



I = I{9), Q = Q{9), U = 0. 



(123) 



for the outgoing wave. This is just the situation with a k-mode of density perturbations at the last 
scattering. As explained in Section 2, The k-mode of density perturbation is azimuthal symmetric 
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about the k axis. At the last scattering, the incident light is unpolarized. Therefore, in Eg. (1151) we 
only need / and Q for a k mode of the density perturbations [2Zl EH [29] , 

Since Q only depends on 6, so that d'^Q = d'^Q = 9^[(1 — ^^)Q{^)], resulting 

E = -dl[{l - ^^')Q{^^)i (124) 

5 = 0, (125) 

i.e., the scalar metric perturbations generate no polarization of magnetic type Another more 

geometric way to see why B = Q and 7^ is to use the definitions in Eg. (11141) . Since U = for 
a k-mode of density perturbation, the polarization matrix in Eg. (II 130 reduces to 

P!r = \i'i _nl..V (126) 



2 V -g sin^ 6 

and, a direct calculation yields 



E = -Q,ee - (127) 
sm d 

2 2cos^ 

sm sm 9 

This tell us that the magnetic type of polarization B essentially involves the derivative of Q with 
respect to (p, and is a measure of asymmetry of polarization field under the rotation about the k 
axis. Since Q is independent of 0, one has i? = 0. 

It is interesting to compare with the case GW, where the rotational symmetry is lost for the k 
mode of GW, and the outgoing light after Thomson scattering would be a general linear polarized 
one, with all three Stokes parameters / = 1(9, (p), Q = Q(9,(f)), and U = U{9,(f)) ^ 0, depending 
on 9 as well as 1211 ES], and resulting \n E ^ and B ^ 0. This distinguished feature of 
non-vanishing magnetic type of polarization of CMB can be served as a possible channel to detect 
gravitational waves. 

The multipole moments af^ of temperature anisotropies and af^ of the electric type of polar- 
ization are given by 

/•I 



aUk) = 2nJ ^dfiYCMHr,^^), (129) 
{l-2)\- 



(/ + 2)! 



df,Y;ME{f^). (130) 



1 



Both af and af are observables on the sky. Since / and E are now functions of 9 only, one can set 
the magnetic index m = in the above expressions and uses the replacements Yio{9) = \J ^^Pi{n) 
and af^ aj^ and af^ ag. 
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Firstly, we calculate the multipole moments at the present time tq. From Eg. (11290 . using 
Eg. (1161) and Eg.(|22]), one has 



al{k) = 27r7 



21 + 1 

A-n 



-1 



27r7 



21 + 1 
47r 



/ 



(i-ff odHi , 1 , , , 



Making use of the relation 



the above expression of is reduced to 



alik) = 7(-2)V47r(2/ + l)ar(fc), 



where 



TO 



dr dr d('^ 

with the variable C = k(To — r). 

Next, we calculate the multipole moments at the present time tq 

^ 121 + 1 '•^ 



dC 



(131) 



(132) 



jKC), (133) 



afo(/c) = 2n 
By Eg. (HZ]) and Eg. ([1211), one has 



(^-2)! 
L(^ + 2)!, 



An 



di2Pi{fi)E{fi). 



Substituting the expression of Eg. (!23l) into Eg. (11351) yields 

3 



(134) 



(135) 



E(/x) 



-7 



4 



dry(r)X2(r)9;l[(l-/.^)^e 



7 



Q /"TO 

^ dTV{T)UT){i+d'^nce-<^). 



(136) 



Substituting Eg. (11361) into Eg. (1 1341) and using Eg. (11310 . one has 

■{1-2) 



af^{k) = 727r 



2^ + 13 

An 4 



{l + 2)\ 



dr / d^^P,{^^)V{r)X2{r){l + ^lf{ee-'^n 



7(-i)V4vr(2/ + 1)- 



(^ + 2)! ' 
(Z-2)!, 



2 /"TO 



rfry(r)X2(r)(l + a,2)^(C^jKC)). (137) 



Using the relation for the spherical-Bessel functions 



j;(x)+2^ + ii-i(^b,w = o 
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to replace j/'(C), the term (1 + d^y{Cji{C)) in Eq. (fT37l) . one obtains 



:i38l 



where 



afik) 







[(^-2)!j 


7 

^0 



(139) 



One can check that of and af in Eqs.f llSSD and (I139p are essentially ak and projected on the 
basis Pz(/i), respectively: 



1 



(140) 
(141) 



The main result of the Appendix is the expressions of af in Eq. fllSSD and af in Eq. (11390 . which 
have been used in Section 5. 
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Figure 1: The visibility function V(t) for the decoupling. The solid lines are given by the analytic formulae 
for different fif, from Ref . |40j . The dash line is the fitting by two pieces of half Gaussian functions as in Eg. (|33|) . 
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Figure 2: The optical depth function k{t) (solid) [40] can approximated by a decreasing exponential function 
cx e"^ I ^'^'i (dots) around the recombination. 
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Figure 3: The multipole moment af has four terms in Ea. (l78|) . which are schematicahy plotted for a 
comparison, ag is dominant at / ~ 200, ai is dominant at I < 100, the ISW is flat and low, and ATdM{Td) 
term is the lowest with two small bumps. 




Figure 4: The perturbation modes /^(t) and 77(t) continuously joined at T2, respectively. 
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Figure 5: The perturbation modes H{t) and Hi{t) continuously joined at T2, respectively. 
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Figure 6: The analytical spectra (red line) TqI{1 + 1)C(^^ /2n compared with the numerical result (dash 
lines) of CAMB [H] and the observed (square dots) WMAP5 [7]. 
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Figure 7: The profile of l{l 
Eq.dSa 

(red dotted line). 



In particular, the bump locations of l{l + l)Cj^^ is determined by that of the scalar perturbations 



(black line) is determined by its source DE{k)M{Td)/ (k{To — Td))^ in 
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Figure 8: It is the time derivative q;o(t) that gives rise to the first two bumps in 
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Figure 9: The spectra C^-^ depend on the prhnordial power index Ug of the scalar metric perturbations. 
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Figure 10: Cf^-^ depend on the running index as of the scalar metric perturbations. 
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Figure 11: 



depend on the baryon fraction fif, 











IN 200 - 










cvj 8000 - 
^ 7000 - 








y: 










^ 6000 - 


















5000 - 








^ 100- 








■i\ 
■M 








/ ^ ■ ■ ^ \ 

/■ ■ I V 


■ o" 






jj 


■I I - 


^ 3000 - 






.'.-A A . 

\ V 


± °~ 








2000- 








H 

-100 - 






^"■'^ A-'l 








\- 












1 









1 




40 














■ — ' 30 








solid: Q^,=0.1 
dash: Q^=0.045 
dots: o =0.01 

b 




K 25 








CNJ 








m 20 

o 

-~- 15 






.■ iJ y. 








...(/ i 








.'// I 
















Figure 12: As the analytic expression tehs, a longer recombination process (greater Ar^) yields a higher 
amplitude of polarization Cf^, and brings more small-scale damping. 
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Figure 14: A late recombination time (larger Td) shifts the peaks and troughs of polarization to larger angular 
scales. 
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Figure 15: The ratio of the analytic spectra to numerical spectra. Left: Cf^ [a) / Cf^ {n) . Right: 
Cf^/(a)Cf^(n). 
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